Spin and charge order in Mott insulators 
and rf-wave superconductors 
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We argue that aspects of the anomalous, low temperature, spin and charge dynamics of the high 
temperature superconductors can be understood by studying the corresponding physics of undoped 
Mott insulators. Such insulators display a quantum transition from a magnetically ordered Neel 
state to a confining paramagnet with a spin gap; the latter state has bond-centered charge order, a 
low energy S = 1 spin exciton, confinement of 5* = 1/2 spinous, and a free S — 1/2 moment near 
non-magnetic impurities. We discuss how these characteristics, and the quantum phase transitions, 
evolve upon doping the insulator into a d-wave superconductor. This theoretical framework was 
used to make a number of predictions for STM measurements and for the phase diagram of the 
doped Mott insulator in an applied magnetic field. 
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I. INTRODUCTION 

Although there is general agreement that the ground 
state of the high temperature superconductors is a con- 
ventional BCS (i-wave superconductor, a number of low 
temperature properties do not fit easily into this frame- 
work, especially in the undcrdoped regime. Among 
these are (a) the presence of a 5* = 1/2 moment near 
non-magnetic Ztl or Li impurties in the undcrdoped 
superconductorEl^cl, (6) low energy collective spin excita- 
tions {a S = 1 spin exciton) at incommensurate wavevec- 
tors near (7r,7r), and (c) instabilities to various co- 
existing spin and charge density wave states ( "stripes" ) in 
certain materials. While it is certainly possible to devise 
microscopic models, and a corresponding Hartree-Fock 
treatment, to generate any of these physical properties 
in the superconducting state, such a procedure is some- 
what ad hoc and has limited predictive power. 

Following an early suggestion of AndersonQ, we explore 
here the idea that these anomalous properties are related 
to fundamental characteristics of the undoped Mott insu- 
lator, and argyue that this leads to a deeper understanding 
of the physicsd. The strategy can be summarized by re- 
ferring to the phase diagram in Fig I. The first step m 
our theoretical program is an understanding of param- 
agnetic Mott insulators which are 'near' the Neel state 
i.e. those obtained across a continuous quantum transi- 
tion at which Neel order is destroyed (this corresponds 
to moving vertically along the 6 = line in Fig |l]) . There 
has been much work on this question in the last decade, 
and we refer the reader to a recent paper by the authorcl 
for a discussion of recent progress along with a review of 
earlier piifffk- The main result we focus on here is the pre- 
dictionoD that the appropriate paramagnetic Mott insu- 



lator has bond-centered charge order and confinement of 
spinous. Fig || contains a sketch of a likely bond-ordered 
state, along with a simple physical argument for its se- 
lection. A correspondingly simple argument for the con- 
finement of 5* = 1/2 spinons in such a state is sketched 
in Fig 1^. The elementary excitation appears naturally 
as a 5 = 1 exciton indicated in Fig |[ (This exciton 
will appear as a 'resonance peak' in the neutron scat- 
tering cross-section, and the first prediction that such a 
peak would appear after the loss of magnetic order in 
doped antiferromagnets (with or without spin-Peierls or- 
der) was made in Ref. 14 ) Another crucial propert}2-nf 
the confining paramagnetic insulator is the predictior£3'cl 
of a free 5 = 1/2 moment near each non-magnetic impu- 
rity: a simple argument for this appears in Fig ^. The 
reader should already notice a strong similarity between 
the properties of the paramagnetic Mott insulating state 
we have just discussed and the properties (a)-(c) of the 
undcrdoped superconductor. We will amplify on this con- 
nection in the following sections, which contain a review 
of recent work, accompanied by a discussion of some rel- 
evant experiments. 



II. NON-MAGNETIC IMPURITIES 

We have already mentioned the recent observation of 
S" = 1/2 moments near non-magnetic impurities at low 
T in the undcrdoped. high temperature superconductors 
by the Orsay groupaB. Specifically, they used NMR mea- 
surements to show that the local susceptiblity near a non- 
magnetic Li impurity has a 1/T divergence as T — > 0. It 
is expected that in this regime the impurity susceptibil- 
ity Ximp(T), formally defined as the excess contribution 
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FIG. 1. Schematic T — phase diagram (adapted from 
Ref. ^) for the high temperature superconductors as a func- 
tion of a ratio of the near neighbor exchange interactions 
and the hole concentration, 5; e.g. the vertical axis could 
be J2/J1, the ratio of the first to second neighbor exchange. 
The shaded region has charge order i.e. at least one of the 
order parameters, cf>x,y, in (EI) is non-zero. The hatched re- 
gion has broken spin-rotation symmetry with (S) 7^ 0, and 
at least one of the order parameters $x,ya in (^) is non-zero; 
(S) — elsewhere. The unfrustrated, insulating antiferro- 
magnet with long-range Neel order is indicated by the filled 
circle. At 5 = 0, there is an onset of charge order above the 
point A, while spin-rotation invariance is restored above B 
(for a review of the S = physics, see Ref. The mag- 
netic ordering quantum transition for 5 > along the line 
BC is expected to be in the same universality class as that 
at t he po int B at 5 = 0, and is discussed as case (I) in Sec- 
The magnetic ordering transition along the line CD 
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difi'ers slightly and is discussed as case (II) in Section 
Neither of these theories apply to the multicritical point C, 
in the, vicinity of which the transition is expected to be first 
orderu. It is not know whether doping the unfrustrated an- 
tiferromagnet (denoted by the filled circle) leads to a system 
which crosses line BC or CD. The nature of the charge orders 
as determined by the computations of Ref. are indicated 
at the top of the figure; numerous transitions, within the gray 
shaded region, in the nature of the charge ordering are not 
shown. The ground state at very low non-zero doping is an 
insulating Wigner crystal and there is subsequently a insula- 
tor-to-superconductor transition; superconductivity is present 
over the bulk of the S > region. The central idea behind 
our approach is that many essential aspects the spin excita- 
tion spectrum of the insulating, paramagnetic region {S — 0, 
(5') = 0) lead to a simple and natural description of the anal- 
ogous properties of the d-wave superconductor. 
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FIG. 2. Representation of antiferromagnetic spin corre- 
lations in the paramagnetic iMott insulator in the language 
of the quantum dimer modelij. Each picture is a snapshots 
of nearest-neighbor singlet valence bonds, represented by el- 
lipses as indicated in (e). The states (a) and (b) (and simi- 
larly, the states (c) and (d)) are connected to each other by 
resonance around a plaquette. The state in (a) has the max- 
imal number of plaquettes around which such resoaance can 
occur, and this lowers its energy. It has been arguedu that this 
effect leads to a broken lattice symmetry in two dimensions, 
and the resulting ground state has bond-centered charge or- 
der. A plausible candidate for such a state is the one with 
spin-Peierls order — this state breaks the square lattice sym- 
metry in the same pattern as the state in (a). The presence 
of this broken symmetry will also modulate the bond charge 
density (i.e. the charge density on an orbital which is a linear 
combination of the orbitals on the two sites at the ends of the 
link) with period 2. 
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FIG. 3. (a) Cartoon picture of the bosonic 5=1 exciton of 
the insulating paramagnet-a similar picture is also expected 
to apply in the superconductor, (b) Fission of the S = 1 exci- 
tation into two S — 1/2 spinous. Notice that the spinous are 
connected by a "string" of valence bonds (denoted by dashed 
lines) which are not able to resonate with their environment: 
this string costs a certain energy per unit length and leads to 
the confinement of spinous. 
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FIG. 4. Confinement ofn a S = 1/2 moment near a 
non-magnetic Zn impuritycS. We consider two Zn impuri- 
ties and imagine moving them infinitely far apart. In the 
trial state in (a), all Cu spins are paired into singlet bonds; 
while this is preferable when the Zn impurities are not too 
far, there is a string of valence bonds connecting them which 
are not able to resonate, and this will eventually lead to the 
preference of the structure in (b). In this case there is no 
string defect in the spin-Peierls order, but each impurity has 
a "free" moment in its vicinity. There is of course a very weak 
residual coupling between these moments even when they are 
far apart, and so the lowest energy states will be a singlet 
and triplet separated by a gap which decreases exponentially 
with the separation between the impurities. This should be 
contrasted with the case when moments do not form around 
each impurity (as may be expected in a deconfined paramag- 
net) when the gap above the ground state saturates at a fixed 
finite value as the impurities move far apart. 
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to the uniform magnetic susceptibility due to a single 
localized deformation (an "impurity" ) obeys 

Xi.p(T^O)^^ (1) 

with C = 1/4. The form (|l|) is very familiar as the 
Curie high temperature limit in a large variety of sys- 
tems. However, we are interested here in the opposite 
T — i- limit, and then can arise only under special 
circumstances. (Note also that in any real experiment, 
there is a finite density of impurities, and even under the 
special circumstances (^ will only hold down to some 
very low temperature below which interactions between 
the impurities have to be accounted for.) We now con- 
sider a wide class of correlated electron ground states and 
list those that can obey (|l|) for an arbitrary localized im- 
purity. 

(i) Fermi liquid: A local spin degree of freedom will obey 
(|l|) only if its exchange interaction with the fermionic 
quasiparticles is ferromagnetic: then we must have C = 
S{S + l)/3 with S integer or half-integer. If the exchange 
is antiferromagnetic (as is surely the case in the cuprates) 
then Kondo screening removes the moment as T and 
Ximp saturates at a finite value. 

(ii) Insulator with Neel order. The magnetic order with 
(S) = picks a preferred direction in spin space, and the 
impurity spin is not free to rotate. The magnetic sus- 
ceptibility is finite only in the direction transverse to the 
Neel order, and the pinning of the spin rotation by the 
magnetic order leads to a finite Ximp fts,T — > 0. So there 
is no S* = 1 /2 moment in a Neel statdlS. 

(iii) Spin gap insulator: As we have shown in Fig U, it is 
easy to find configurations (even with non-magnetic im- 
purities) which have C = S{S + l)/3 with S intcgcr/half- 
integer. 

[Wl— Magnetic quantum critical point: We recently stud- 
iedtj the critical point between (ii) and (iii) in two dimen- 
sions. We found there that (|l|) continues to be obeyed 
but C is an irrational number — thus, remarkably, the im- 
purity is associated with an irrational spin, 
(v) d-wave superconductor. This case similar to the Fermi 
liquid in that it is possible to have Kondo screening of 
the moment as T — > 0. However th^-conditions for Kondo 
screening are much more stringentEII because of the lower 
density of fermionic states at low energy: it is only 
present if the antiferromagnetic exchange and the degree 
of particle-hole asymmetry exceed some finite value. In 
the free moment phase C — S{S + l)/3 with S integer or 
half-integer, while at the impurity quantum critical point 
between the free momeat and Kondo screened phases C 
is an irrational numbered. 

The NMR observations show that there is no Kondo 
screening of the moment at low doping, and that there 
is an onset of a non-zero Kondo temperature (Tk) at a 
critical doping beyond which the Tk rises rapidly. These 
results can be very naturally explained in a model of a 
5=1/2 moment interacting with a di-waAie superconduc- 
tor, as has been shown in recent workll5t3. The separate 



question of how the underdoped superconductor has a 
moment near a non-magnetic impurity in the first place 
is answered for "free" by the theoretical framework of 
Fig 0. Moving along the insulating vertical line at (5 = 
in Fig |l| we proceed from case (ii) to case (iii) via (iv) : a 
precise theory of the fornmtion of the moment along this 
route has been presentedlij. Then we dope the param- 
agnetic charge-ordered insulator to reach the supercon- 
ductor with (S) — 0; along this path the question is not 
one of the formation of the moment, but instead of its 
quenching by the fermionic quasiparticles; the latter is- 
sue is precisely that addressed by the Kondo model. Full 
translational symmetry must be restored at some bulk 
quantum critical point along this path (to be discussed 
in Section III), but this point will not, in general, coin- 



cide with the impurity quantum critical point at which 
Tk first becomes nonzero. 

Next weJjurn to STM measurements of non-magnetic 
impuritiesEil: these show a sharp low bias peak in the 
tunneling conductance at very low T in optimally doped 
samples. It was argued in Ref. 18 that this peak is related 
to a Kondo resonance in the rf-wave superconductor and 
so its bias and width are determined by Tk ■ An immedi- 
ate prediction is that there should be dramatic change in 
the bias dependence and spatial distribution of the tun- 
neling conductance as the T is increased through Tk- In 
particular, the states near the impurity cross over from a 
coherent Kondo resonance to ones determined by poten- 
tial scattering off the chemically distinct orbitals on the 
Zn site. 

A more subtle issue, which is far from resolved, is that 
of the spatial dependence of the low-bias peak in the low 
T STM. In Ref. ^the author and collaborators proposed 
that the unusual spatial dependence may be explained 
by Kondo scattering from a spatially diffuse S" = 1/2 mo- 
ment on the four nearest neighbor sites of the Zn moment. 
However, this computation did not account .fep-the strong 
potential scattering from the central Zn siteE3Ej, and also 
the possible 'blocking' or 'filter' effect of the interveai 
BiO layer between the STMJjip and the Cu02 layeiuJ' 
The work of Zhu and Tinepj, and our discussions above 



and in Section [II instead suggest the following model 
for future study: over the time scale of an electron tun- 
neling event it is possible that local "stripe" correlations 
preferentially pick out a specific nearest neighbor site of 
the Zn irnpurity upon which the moment resides, just 
as in Fig ^ . A Kondo model with the moment localized 
on this site, along with a strong repulsive potential on 
the Zn site, should then be solved for its local density of 
states. The actual STM measurements, of course, con- 
stitute a long time average of the tunneling current and 
this is obtained from the theory by an incoherent average 
of the results of the four possible positions of the local 
moment. Note that it is the relatively long-lived stripe 
order which demands an incoherent average, in contrast 
to the coherent average used in Ref. |l^ which neglected 
stripe correlations. 
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III. SPIN AND CHARGE ORDERING 
QUANTUM TRANSITIONS 

We will only consider spin and charge order in an ideal- 
ized doped antiferromagnet, where any random disorder 
potential generated by the dopant ions is neglected. In 
the presence of disorder, the spin order (and oancomitant 
charge order) may appear as spin glass ordeitZI. Related 
ideas have been discussed by ZaanennS. 



A. Charge order 

We begin by reviewing simple Landau-theory-like con- 
siderations ta characterize charge order in doped anti- 
ferromagnetsQ. We consider charge density wave fluctu- 
ations in the (1,0) and (0,1) directions for a lattice in 
which the Cu ions occupy the sites at integer values of 
X, and the O ions occupy half-integer x. Then we can 
describe the spacetime (r ~ ix,y)) modulations in the 
charge density, Sp{r,T) by 



Sp{r, t) 



4>x{r,T) + c.c. + X 



V 



(2) 



where Kc is the ordering wavevector and (l)x,y are complex 
order parameters. For commensurate Kc = 2-Kmjn (m, 
n, integers) the coupling to the underlying lattice will 
allow a term in the action 



Ac / (frdT(t>2{r,T) + c.c. -f ; 



y- 



(3) 



This term pins the phase of (j)x,y either at arg[(/)j; ,^] ~ 
for Ac < 0, or at arg[(/)^_j^] — tt/ti, for Ac > (changes 
of arg[0a;,y] by integer multiples of 27r/n are not signif- 
icant as they correspond simply to a translation in (|^) 
by integer lattice spacing); the first case corresponds to 
a site- centered charge density wave, while the second is 
bond- centered. The state in Fig|2^ has m = 1, n = 2 and 
is bond-centered; note that for this case only the charge 
density on all the Cu sites in unmodified — the modula- 
tion is apparent if we examine the charge density on the 
O sites, or in off-site correlations of the Gu-orbitals which 
take the form of a px density- wave statecj. 

To choose between the different possibilities for (||) 
we have to turn to more microscopic considerations. In 
the regiQn with (S) = in Fig ^ a controlled large N 
theoryuti was used to follow the evolution of the 6 = 
bond-charge order in Fig|^a to non-zero 6: we found only 
bond-centered waves with n even. The value of n jumped 
to a large number at small d and then decreased mono- 
tonically through a series of even integer plateaus with 
increasing 5 (see the top of Fig |l]) ; in particular there was 
large plateaus at n = 4 (and occasionally a small plateau 
at n = 2) before translational symmetry was eventually 
restored to an isotropic d-wave superconductor. (Spe- 
cific theories of hole doping and icrmionic excitations in 
states with n^-f= 2 charge orderM, and also with "pla- 
quette" orderE3, have been presented recently.) We have 



proposed that such bond-centered charge order with even 
n is characteristic of the high temperature superconduc- 
tors in the region with no spin order i.e. (S) = 0; in the 
regime where the ground state is an isotropic d-wave su- 
perconductor (see Fig |]) strong low energy fluctuations 
of such bond charge order can be expected. Some exper- 
imental support for this picture has emerged in recent 
measurements of the phonon spectraOO. 

The nature of the quantum theory describing the even- 
tual restoration of translational invariance has also been 
discussedcl: the critical fixed point between a d-wave 
superconductor and a d-wave superconductor with co- 
existing charge order is described by a relativistic field 
theory with dynamic exponent z = 1 which obeys hyper- 
scaling and exhibits w/T scaling. This shpuld be con- 
trasted with the work of the Rome groupE3 who have 
considered charge ordering transitions in Fermi liquids: 
such fixed points are very different and do not satisfy any 
of the characteristics just mentioned. 



B. Spin order 

A related analysis can be developed for spin excitations 
in the doped Mott insulator. Here it is convenient to be- 
gin in the paramagnetic Mott insulator at S = which, 
as shown in Fig ||, has a sharp S' = 1 exciton above a 
spin gap. Upon moving to the doped superconductor, it 
is easy to see that this exciton retains its integrity for 
a finite range of S. The superconductor does have an 
entirely new class of low energy spin excitations which 
are not present in the insulator: these are the fermionic 
5=1/2 Bogoliubov quasiparticles. However, these low 
energy quasiparticles reside only in certain sectors of the 
Brillouin zone and momentum conservation constraints 
can prevent the S = 1 exciton from decaying into two 
5=1/2 quasiparticles. Explicit calculations exhibit- 
ing this feature appear in Ref. |ll|. So the main effect 
of the non-zero 6 will be renormalize the dispersion of 
the exciton: in particular, it is likely that the exciton 
will eventually move its minimum away from (tt, tt) to 
a general wavevector (tt, vr) + Kg. As in the discussion 
of the charge density, we assume that Kg is polarized 
along (1,0) or (0,1). In the presence of such excitonic 
fluctuations, the spin density on the Cu sites is given by: 

5c(r, r) = (-l)-+^e^^=^$,c(r, r) -t- c.c. +x^y (4) 

where a = x,y, z are the spin directions, and ^xa^ ^ya 
are complex fleld operators for the 5 = 1 excitons. We 
now discuss various terms in the effective action for i 
^ya. For the reasons just discussed, we will neglect their 
couplings to the fermionic quasiparticles, but it is not 
difficult to include such a coiipjing if momentum conser- 
vation happens to permit itBE3. First, along the lines of 
Refs. |,|l|,|| we have the usual gradient and time deriva- 
tive terms present in any quantum antiferromagnet: 
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S<h = d rdr 



Ul 



|2\2 I /irTv |2\2 



(5) 



Second, the underlying lattice may pin Ks to the com- 
mensurate values Ks = or Kg = 27r/(2p) withp integer. 
For the special case Kg — 0, we see from (^) that '^x.y 
have the same spin modulation and so only one of them 
need be considered; further they should also be real, and 
so 5$ reduces to the familiar 0(3) ip'^ field theory in 2+1 
dimensionsEa. For Kg — 7r/p, as in (^, we have a term 
in the action 



Sl = Xs d^rdT{<^l^f + C.C. 



(6) 



which again pins the phase of '^xa and to a discretje 
set of allowed values. Finally, as noted by Zachar et al\!\ 
there is also a coupling between the spin and charge de- 
grees of freedom 



5< 



$0 



w 



d^rdT(t>l<^l 



(7) 



which is operative only if Kc — 2Ks and then pins the 
long-range spin and charge order to each other. The or- 
dering observed in LNSCO is of this type with to = 1, 
n = 4, and p = i. Spin correlations pinned at p = 4 
are also universally observed over a wide doping range 
in the cuprates: this pinning is difficult to understand 
from the semiclassical theory of stripe formatioiO, and 
we propose that it is linked via (|^) to singlet bond charge 
order correlations for which a large plateau at n — A was 
found in Rcf. |9[ 

The above framework can now be used to study the 
transition into the magnetically ordered region with 
(S) 9^ in Fig ||. This transition is driven by decreasing 
the value of s in (|^) and corresponds to a condensation of 
the S = 1 exciton. The nature of this transition depends 
crucially on the role of the charge order, and two distinct 
cases can be separated. 

(I) In the first case, long-range charge order remains an 
innocent spectator to the spin ordering transition. If 
the spin ordering is at a wavevector Kg — Kc/2, then 
the charge ordering, via the coupling (Q), will select ^xa 
over ^ya (say) , and also pin its phase so that ^xa is real. 
The theory reduces to that of a single real S-componect 
vector, and is again the familiar 0(3) (p^ field theorya. 
This situation is the case along the BC portion of the 
spin-ordering phase boundary in Fig |l|. 

(II) In the second case, the spin-ordering fluctuations are 
paramount, and charge order appears only in the region 
with (S) 7^ (assuming the spin correlations are not spi- 
ral); this is the case along the CD portion of the phase 



boundary in Fig 0. The critical theory is now given by 
5$ -|- Sl and the spectator charge order is determined 
simply by (^ as (j)x ^xan and similarly for 4>y. 

Finally, we note the nature of the dynamic spin fluctu- 
ation spectrum near the spin-ordering transition in the 
region with (S) = 0. Here, an adequate picture emerges 
if we simply treat (^ in the Gaussian approximation, 
with s > 0. A combination o£-fi) and ^ then implies 
the dynamic spin susceptiblityO: 

^^^'""^ clikx - Gl.,f -h cliky - Glyf + . - (c. + ^r^f 



cf{kx 


- G-'-xY 


+ ^2(^1/ 


- G%Y - 
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- s ~ {uj + iviY 


C\{ky 


~ ^+y) 


+ C2 (kx 


-G\x?^ 
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- s — {oj + iviY 


C\{ky 


~G\f 


+ c^(kx 


-Glx?- 


- s — + iviY 



(8) 



where = (tt ± Xs,7r) and = (tt, tt ± Ks), and 77 
is a positive infinitesimal. A related, spectrum has been 
discussed recently by Batista et aZHj but on the ordered 
side of the transition: they have argued that the super- 
position of the terms in (|^) leads to a maximum in the 
spectral density at the commensurate point k = (tt, tt) 
and this is responsible for the "resonance peak" in neu- 
tron scattering. 



IV. EFFECT OF AN APPLIED MAGNETIC 
FIELD 

In this final section we briefly review recent work0 on 
the effect of an applied magnetic field, oriented perpen- 
dicular to the Cu02 layers, on the above spin and charge 
ordering transitions within a superconducting ground 
state. 

So far, we have found that there was little mate- 
rial difference between the spin/charge ordering transi- 
tions of a superconductor and those of an insulator: the 
superconductor has additional low energy quasiparticle 
excitations, but these often decouple from the critical 
spin/charge degrees of freedom because of constraints 
from momentum conservation. The situation changes 
dramatically in the presence of an external field. Briefly 
stated, the "background" superconducting order has an 
infinite diamagnetic susceptiblity and this feeds into an 
anomalously strong response of the spin/charge order, 
much stronger than that would be found for the corre- 
sponding quantum transitions in an insulator. The effect 
is quite simple: in the presence of an applied fleld, iJ, 
the type II superconductor reacts by setting up superflow 
currents which screen the field. A standard calculation 
using the standard Ginzburg Landau free energy shows 
that 
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mean kinetic energy density of superflow • 



H 


In 


Hc2 




Hc2 




H 



(9) 



This shift leads to concomitant-phange in the effective 
action for the spin/charge orderEZi via the simple Landau 
theory coupling (emphasized in Ref. pSl) 



5, 



K / d2rdT|i/;|2|$^„|2+a;^y 



(10) 



(and similarly for (f>x,y) to the superconducting order 
parameter tp. Crudely speaking, this correction is ac- 
counted for an i/-dfiBendent shift in the value of s in (|^) 
which has the formEZI 



H 


In 


Hc2 




Hc2 




H 



(11) 



This simple expression now has a number of strong conse- 
quences for the position of the phase boundary for the on- 
set of long-range spin order, and for the nature of the spin 
excitation spectrum on either side of the phase boundary. 
On the side with long-range spin order, ( pT| ) implies that 
the Bragg-eJastic scattering intensity should increase (for 
K > 0) a£3 ^ (|iJ|/i/c2) ln(i?c2/|-ffW: this prediction is 
consistent with recent experimentsE3. On the side with 
no spin order, the dynamic susceptibility will be given 
by (^ but with s replaced by s{H): hence the minimum 
excitonic energy, yJs{H), will decrease with inciieasing 
H — this trend is also consistent with experimentald. 

All of the above effects considered the uniform, spa- 
tially averaged, consequence of the applied field. How- 
ever, as is well known, the superflow kinetic energy in (|^) 
is inhomogenously distributed in the form of a vortex lat- 
tice, and this has interesting observable consequences for 
the spin/charge order. In the spin^rdered phase, there 
are weak satellite Bragg peakgid'EZi around G^'^, sepa- 



m around G^'^ 
rated from G^'^ by reciprocal lattice vectors of the vortex 

lattice. In the region without spin order {{S) = 0) which 
was considered in Ref. |3^, the wavefunction, '^x,ya{r), of 
the lowest energy exciton (which is an excited state at 
energy s{H)) will be peaked at the vortex core, but 
will exte nd w ill outside the vortex core to a distance 
^ ci,2/\/ s{H) (up to the natural maximum of the vortex 



lattice spacing), where the velocities ci_2 are expected to 
be of order the spin-wave velocity in the antiferromag- 
netic state at (5 = 0. Given the large size of the exciton 
wavefunction, small local disorder may well enhance 
over (or vice versa) in a given set of vortices. Associ- 
ated with this excitonic excited state is a spin modulation 
specified by ^x,yair) and (^), and a charge modulation 
given by <j)x,y(.r) ~ ^l^y^ir) and (|). 
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